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Molecular  Seale  Theoretical  Studies  of  Energy  Deposition  and  Redistribution  in  Crystalline 
i  High  Explosives  to  Stimulate  Enhanced  Detectable  Signatures 

Thomas  D.  Sewell 

Department  of  Chemistry;  University  of  Missouri-Columhia;  Columbia,  MO  6521 1-7600 

In  this  project  we  have  developed  two  molecular  dynamics-based  methods  for  simulating  the 
THz-region  infrared  (1R)  spectrum  of  molecular  crystals.  Both  methods  were  designed  for  use  at 
finite  temperature  and  stress.  The  first  method,  referred  to  as  the  mode-relaxation  approach ,  requires 
(1)  the  normal  mode  (/.<?.,  phonon)  eigenvectors  and  frequencies  of  the  cry  stal,  (2)  the  associated 
squared  transition  dipoles  p2  for  all  THz-region  phonons  with  zero  wave  vector  (k  =  0),  and  (3)  the 
spectral  line  width  y  for  each  k  =  0  THz-region  phonon  with  nonzero  p2  (i.e.,  the  lR-active  THz- 
region  modes).  These  quantities  are  the  ones  required  to  specify  the  Lorentzian  absorption  profile  for 
an  individual  spectral  line.  The  eigenvectors,  frequencies,  and  transition  dipoles  were  obtained 
directly  from  normal  mode  analysis  of  the  crystal,  performed  for  crystal  lattice  parameters 
corresponding  to  a  specified  temperature  and  stress  state.  Line  widths  were  determined  separately  for 
each  THz  mode  by  selectively  depositing,  in  an  otherwise  thermalized  crystal,  a  fixed  amount  of 
energy  in  a  given  normal  mode  and  monitoring  the  time  dependence  of  energy  transfer  from  that 
mode.  Relaxation  time  constants  for  each  mode  were  determined  by  fitting  the  corresponding  energy 
decay  time  history,  averaged  over  an  appropriate  ensemble  of  trajectories,  to  an  exponential  function; 
from  which  the  line  width  was  obtained  using  y  =  l/(47tcx).  The  Lorentzian  profiles  for  the  IR-active 
modes  were  summed  to  yield  the  absorption  spectrum.  The  advantages  of  the  mode-relaxation 
approach  are  that  it  yields  (1)  detailed  information  concerning  which  particular  normal  modes 
contribute  to  the  THz  absorption  intensity  and  (2)  time  scales  and  pathways  for  energy  transfer  from 
the  selectively  excited  THz  modes  into  the  remaining  modes  of  the  system. 

In  the  second  method,  referred  to  as  the  DDACF  approach,  the  1R  absorption  intensity  is 
obtained  as  the  Fourier  transform  of  the  dipole-dipole  autocorrelation  function  of  the  crystal 
calculated  from  an  ensemble  of  equilibrium  trajectories  for  a  given  pressure  and  temperature.  In 
contrast  to  the  mode-relaxation  approach,  the  DDACF  approach  does  not  provide  information  about 
particular  mode-mode  energy  transfer  pathways.  It  does,  however,  sample  directly  the  effective 
modes  of  the  system,  as  opposed  to  the  set  of  zeroth-order  normal  modes  obtained  from  normal  mode 
analysis,  at  any  temperature  or  for  any  defect  structure  in  the  crystal. 

Both  methods  were  evaluated  by  predicting  the  THz  spectrum  of  the  secondary  explosive 
pentaerythritol  tetranitrate  (PETN)  for  a  variety  of  thermodynamic  conditions.  Specifically,  using  a 
force  field  known  to  accurately  predict  several  thermodynamic  and  mechanical  properties  of 
crystalline  PETN,  simulations  were  performed  for  defect-free  crystals  at  T  =  298  K  and  hydrostatic 
pressures  of  P  =  0,  1,2,  and  3  GPa.  Simulations  were  also  performed  at  T  =  298  K  and  P  =  0  GPa  for 
crystals  containing  vacancy  defects.  Additional  studies  were  performed  for  defect-free  cry  stals  at  T  = 
298  K  and  P  =  0  GPa  to  assess  the  sensitivity  of  the  predictions  to  modest  modifications  of  the  force 
field  model;  namely,  replacement  of  harmonic  oscillator  covalent  bond-stretching  terms  by  Morse 
oscillators,  or  application  of  geometric  constraints  to  fix  C-H  and  N-0  covalent  bond  distances  at  the 
corresponding  equilibrium  values.  Both  methods  predict  IR  spectra  at  T  =  298  K.  and  P  =  0  GPa  that 
are  in  good  agreement  with  experimental  data,  although  the  predicted  peak  widths  are  noticeably 
narrower  than  the  observed  ones.  There  are  no  high-pressure  data  available  for  direct  comparison. 
However,  the  predicted  THz-region  1R  pressure  coefficients  dco/dP  agree  well  with  Raman  data 
obtained  for  slightly  higher  frequencies.  The  two  principal  effects  of  vacancy  defects  on  the  predicted 
spectra  are  red  shifting  of  the  main  absorption  peaks  and  broad,  increased  absorption  intensity  at  the 
lower  end  of  the  frequency  range. 
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